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SECTION  I 


INTRODUCTION 

The  Multi-Hundred  Watt  Radioisotopic  Thermoelectric  Generator  (MHW-RTG)  is 
used  as  the  power  generation  system  for  the  NASA  Voyager  1977  spacecraft.  The 
use  of  plutonium  fuel  In  these  RTGs  necessitated  an  extensive  nuclear  safety 
evaluation  and  risk  assessment  of  the  Voyager  launches.  This  safety  evaluation 
consisted  primarily  of  a review  of  the  RTG  contractor's  Final  Safety  Analysis 
Report  (FSAR)  (ref.  1)  and  preparation  of  a Safety  Evaluation  Report  (SER)  by 
the  Interagency  Nuclear  Safety  Review  Panel  (INSRP).  One  area  of  concern  which 
received  a great  deal  of  attention  during  the  safety  review  was  the  possibility 
of  an  accidental  reentry  of  the  spacecraft  into  the  earth's  atmosphere. 

The  total  reentry  analysis  was  a complex  task,  involving  several  organiza- 
tions. The  General  Oynamics/Convalr  Division  calculated  the  initial  vehicle 
velocity  and  entry  angle  (V-^)  conditions  (ref.  2).  From  this  V-Y  map,  four 
reentry  conditions  were  selected  for  detailed  analysis:  shallow,  maximum  stress, 
orbital  decay,  and  powered  reentry.  The  Jet  Propulsion  Laboratory  studied 
spacecraft  breakup  and  determined  the  conditions  at  the  time  of  RTG  release  for 
each  of  these  reentry  conditions  (ref.  1,  Appendix  B).  The  General  Electric 
Space  Division  (the  RTG  contractor)  analyzed  the  RTG  and  Heat  Source  Assembly 
(HSA)  thermal  responses,  the  release  of  the  HSA,  and  the  aeroshell  and  ablation 
sleeve  stresses  (ref.  1,  Appendix  E).  It  was  found  that  the  maximum  stresses 
in  the  aeroshell  occurred  for  a V-y  of  36,000  fps  and  -46"  (at  400,000  ft)  and 
at  12.9  s after  reentry  (defined  as  400,000  ft). 

If  the  maximum  allowable  stress  was  exceeded,  the  aeroshell  would  fail 
causing  release  of  the  internal  Fuel  Sphere  Assemblies  (FSAs);  therefore,  a 
thorough  stress  analysis  of  the  reentering  HSA  was  required.  The  preliminary 


1.  Final  Safety  Analysis  Report  for  the  MJS-77  Mission,  Doc.  No.  77SDS4206, 
General  El ectHc  Company , Space  Division,  Ph i 1 adel phia , PA,  January  1977. 

MJS  '77  RTG  Safety  Study  - Phase  II.  Range  Safety  Equipment,  Launch  Pad 
Hazards,  Launch  Vehicle  Failure  Probabilities  and  Reentry  Environment, 

Doc.  No.  CASO/LVP  ?S-004,_  General  bynamics  Convair  Division,  San  Diego,  CA, 
7 July  1975. 
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analysis  Involved  a plane  stress  approximation  for  the  aeroshell  and  ablation 
sleeve  center  sections  using  the  GE  ORTHOSAFE  computer  program.  This  approxi- 
mation neglected  axial  stresses,  so  the  second  technique  modified  the  plane 
stress  analysis  for  three-dimensional  effects  by  applying  a longitudinal  stress 
distribution  to  maintain  a plane  aeroshell  cross  section.  The  final  analysis  was 
a three-dimensional  stress  calculation  using  the  SAP  IV  computer  program.  This 
analysis  used  an  R-0  temperature  profile,  but  assumed  the  temperature  was  constant 
along  the  aeroshell  length.  (For  a detailed  discussion  of  the  stress  analyses, 
see  reference  1,  Appendix  E.) 

The  results  of  these  three  analyses  are  shown  In  Table  1.  The  maximum 
allowable  aeroshell  stress  was  7214  psl,  and  therefore  It  was  not  predicted  to 
fail;  l.e.,  a positive  margin  of  safety  was  calculated. 

Table  1 

RESULTS  OF  THE  AEROSHELL  STRESS  ANALYSIS 

Plane  Stress  3-D  Modification  SAP  IV 

Maximum 

Calculated  6050  6949  6186 

Stress 


The  plane  stress,  three-dimensional  modification,  and  SAP  IV  calculations 
described  above  neglected  the  effect  of  a three-dimensional  temperature  field. 
This  report  Investigates  this  effect,  specifically  on  the  SAP  IV  results.  Since 
the  longitudinal  temperature  variations  are  not  known  at  12.9  s after 
reentry,  the  precise  effect  on  the  stress  magnitudes  and  distributions  cannot  be 
determined;  however,  the  three-dimensional  temperatures  are  known  at  13.3  s after 
reentry.  Combining  the  temperatures  at  these  two  different  times  during  the 
reentry  would,  of  course,  be  meaningless.  But  If  the  temperature  gradients 
along  the  aeroshell  length  are  on  the  same  order  of  magnitude  for  the  two 
slightly  different  times,  then  an  estimate  of  the  change  In  predicted  aeroshell 
stress  at  12.9  s due  to  three-dimensional  temperatures  can  be  calculated  by 
applying  the  13.3  s longitudinal  temperature  gradients.  The  results  of  the 
SAP  IV  calculations  presented  below  are  therefore  hypothetical  axial  temperature 
gradients,  and  can  only  be  used  for  order-of-magnitude  estimates. 
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A brief  description  of  the  Input  mesh  and  boundary  conditions,  material 
properties,  temperature  profiles,  and  external  loadings  will  be  made  (largely 
taken  directly  from  the  FSAR),  followed  by  a presentation  of  the  SAP  IV  results 
and  the  conclusions. 
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SECTION  II 

THREE-DIMENSIONAL  FINITE  ELEMENT  MESH 

A schematic  diagram  of  the  three-dimensional  mesh  used  as  input  to  the 
SAP-IV  program  is  shown  in  figure  1.  This  mesh  is  identical  to  that  used  in  the 
FSAR.  It  contains  1776  elements  and  2475  node  points. 

Boundary  Conditions:  Due  to  symmetry,  only  one-quarter  of  the  aeroshell 
was  modeled.  The  nodes  lying  in  the  X-Y  centerplane  were  restrained  in  the 
Z-direction,  and  the  nodes  lying  In  the  X-Z  centerplane  were  restrained  in  the 
Y-dlrection.  The  third  boundary  condition  was  artificially  introduced  to  simu- 
late the  effect  of  the  endcap. 

Although  the  aeroshell  endcap  was  not  modeled  in  this  analysis,  its  contri- 
bution in  maintaining  the  shape  of  the  end  of  the  aeroshell  was  simulated.  This 
rigidity  effect  was  modeled  by  applying  an  Internal  load  at  the  lock-ring  groove. 
This  load  was  varied  in  the  X-Y  plane  (or  R-0  direction),  and  its  profile  was 
opposite  in  direction  to  that  of  the  opposing  aerodynamic  load  (see  Section  V). 

A scaling  factor  was  applied  to  the  Internal  load  and  this  factor  was  adjusted 
to  restrict  movement.  (Since  a costly  and  time-consuming  Iteration  technique 
was  used,  exact  counterbalancing  of  the  loads  was  impractical.  The  iterations 
were  therefore  terminated  when  the  change  in  aeroshell  stresses  was  less  than 
1 percent.) 
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SECTION  III 
MATERIAL  PROPERTIES 

Temperature  dependent  material  properties  were  used  for  the  ae^oshell  POCO 
graphite  AXF-5Q.  The  elastic  modulus,  Poisson's  ratio,  and  linear  coefficient 
of  expansion  are  shown  in  figures  2,  3,  and  4,  respectively.  These  curves  were 
taken  from  the  FSAR.  The  elastic  modulus  curve  is  an  average  based  on  the  data 
obtained  from  an  earlier  program  (Lincoln  Experimental  Satellites  8 and  9 aero- 
shell). 
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Figure  2.  Initial  Elastic  Modules  vs  Temperature 
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Figure  3.  Poisson’s  Ratio  vs  Temperature 


Figure  4.  Linear  Coefficient  of  Expansion  versus  Temperature 
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SECTION  IV 

TEMPERATURE  PROFILES 

The  two-dimensional  aeroshell  temperature  profile  was  calculated  and 
reported  in  the  FSAR  using  the  THTDA  (Transient  Heat  Transfer  - Version  D - 
Ablation)  code.  This  temperature  distribution  was  calculated  for  R-0  spacial 
dependence,  and  was  assumed  constant  along  the  aeroshell  axis  (or  Z-dlrectlon) . 
The  R-0  temperature  profile  given  for  12.9  s after  reentry  Is  shown  in  Table  2 

Table  2 

AEROSHELL  TEMPERATURE  MATRIX  (Reference  1) 


9,  Degrees 


R,  in. 

0.0 

21.0 

42.5 

66.4 

90.0 

118.6 

151.5 

18.0 

3.280 

2243 

2433 

2380 

2244 

2040 

2021 

2013 

2006 

3.410 

2584 

2663 

2596 

2382 

2105 

2027 

2015 

2010 

3.535 

3109 

3117 

3057 

2681 

2218 

2045 

2026 

2021 

3.595 

3463 

3450 

3442 

2937 

2314 

2064 

2038 

2033 

The  three-dimensional  aeroshell  temperature  distribution  was  not  calculated 
for  the  time  after  reentry  for  maximum  combined  stress  (thermal  plus  aerodynamic), 
which  corresponds  to  a reentry  time  of  12.9  s after  400,000  ft.  The  ax^al 
distribution  is,  however,  available  for  the  time  at  maximum  thermal  stress 
(13.3  s)  from  the  end  section  analysis.*  This  analysis  used  the  THTDA  code  to 
calculate  the  circumferential  temperature  distribution  in  the  aeroshell  endcap. 
Although  It  is  unrealistic  to  blend  the  two-dimensional  profile  at  12.9  s with 
the  axial  distribution  at  13.3  s,  the  effect  a three-dimensional  temperature 
field  has  on  the  aeroshell  stresses  can  be  demonstrated  (however,  only  quantita- 
tively estimated)  by  applying  this  temperature  gradient  in  the  axial  direction. 


♦Letter  from  G.  Drenker  (MJS  Program  Manager),  GE  Space  Division,  to  C.  P.  Melfl, 
Aerospace  Corporation,  El  Segundo,  CA,  10  March  1977. 
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The  axial  temperature  gradients  are  plotted  in  figures  5 through  8.  The 
indicated  gradients  were  applied  to  the  temperature  profiles  in  Table  2 
beginning  at  2.5  in  from  the  end  of  the  aeroshell . From  the  centerplane  to 
2.5  from  the  end,  the  temperature  was  assumed  constant  in  the  axial  (Z)  direc- 
tion. 
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Axial  Temperature  Distribution  at  a Radius 


al  Temperature  Distribution  at  a Radius 
3.520  in  at  13.3  s after  400  kft 


Axial  Temperature  Distribution  at  a Radius 
of  3.370  in  at  13.3  s after  400  kft 


Temperature  Distribution  at  a Radius 
60  in  at  13.3  s after  400  kft 
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SECTION  V 

AEROOYNAMIC  AND  INERTIAL  LOADING 

» 

The  aerodynamic  and  Inertial  loads  at  12.9  s were  applied  to  the  aeroshell 
as  shown  in  figure  9.  These  loads  are  as  described  in  the  FSAR  (ref.  1).  The 
aerodynamic  load  varied  as  a Pcos2o  distribution,  where  P is  the  stagnation 
pressure  (P  * 118  psi  (ref.  1)).  The  inertial  loads,  due  to  Internal  forces  of 
the  fuel  sphere  assemblies,  were  constant  along  the  indicated  arcs.  The  aero- 
dynamic and  Inertial  loads  were  applied  uniformly  along  the  length  of  the  aero- 
shel 1 . 
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SECTION  VI 
RESULTS 

The  mesh,  material  properties,  temperature  fields,  and  external  loads  des- 
cribed above  were  Input  to  the  SAP  IV  computer  program  (refs.  3 and  4).  The 
code  was  run  for  two  cases,  to  calculate  thermal  stresses  only  and  then  to  cal- 
culate thermal  stress  plus  aerodynamic  and  Inertial  loads.  These  two  calculations 
were  made  for  each  of  the  two  temperature  fields,  first  for  the  case  of  a con- 
stant temperature  along  the  aeroshell  length,  and  then  for  the  case  of  an  axial 
temperature  grad i ent . 

The  results  of  these  calculations  are  shown  in  figures  10  through  15.  The 
critical  area,  where  maximum  tensile  stress  occurred,  was  along  the  stagnation 
line  at  the  inside  diameter.  Figures  10  and  11  show  these  stress  profiles. 

Along  the  Inner  diameter,  the  maximum  stress  is  seen  to  occur  at  an  area  6.25 
to  6.50  in  from  the  aeroshell  centerplane,  just  above  the  lock-ring  groove. 

The  stress  through  the  thickness  of  the  aeroshell  at  this  point  is  shown  in 
figures  12  and  13,  and  the  stress  around  the  aeroshell  inner  circumference  is 
shown  in  figures  14  and  15.  The  results  are  summarized  in  Table  3. 


Table  3 

RESULTS  OF  THE  CALCULATIONS,  SHOWING  THE  INCREASE 
IN  MAXIMUM  STRESS  WITH  AN  AXIAL  TEMPERATURE 


Maximum  Combined 

Principle  Stresses 

Percent 

KSI 

Increase 

Constant  Axial  Temperatures 

7066 

10.8 

Axial  Temperature  Gradients 

7828 

3.  Bathe,  K.  J.;  E.  L.  Wilson;  F.  E.  Peterson;  SAP  IV  - A Structural  Analysis 
Program  for  Static  and  Dynamic  Response  of  Linear  Systems,  EERC  73-11 
University  of  California,  Earthquake  Engineering  Research  Center,  Berkeley, 
CA,  June  1973. 

4.  Melfl,  C.  P.;  J.  L.  Johnson;  Viking  and  LES  8/9  Reentry  Nuclear  Safety 
Studies,  Volume  II,  Supplementary  Reentry  Response  Analysis,  aPWL-YR-76-1 64 
Vol  II,  Air  force  Weapons  Laboratory,  Kirtlana  Air  Force  ftase,  NM, 

September  1976. 
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Stagnation  Line  for  the  Constant  Axial  Temperatures 
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Figure  13 


. Stress  Profile  Through  the  Aeroshell  Thickness  at  the  Stagnation 
Line  for  the  Applied  Gradients  to  the  Axial  Temperatures 
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SECTION  VII 

/ 1 CONCLUSIONS 

S' 

This  report  has  demonstrated  the  increase  In  calculated  stresses  when  a 
longitudinal  gradient  is  applied  to  an  R-e  temperature  field.  The  results 
indicate  an  Increase  in  calculated  stress  when  an  axial  temperature  gradient 
is  applied,  and  if  the  magnitudes  of  the  gradients  at  12.9  and  13.3  s after 
* reentry  are  similar,  an  increase  of  approximately  10  percent  is  Indicated. 
Reference  5 describes  an  alternative,  probabilistic  (error  analysis)  approach 
which  may  be  useful  in  the  risk  assessment  for  accident  analyses  which  indicate 
decreasing  margins  of  safety. 


Anderson,  0.  C.,  et  al..  Voyager  - Analysis  and  Risk  Assessment 
77-161,  Air  Force  Weapons  Laboratory,  Kirtland  Air  Force  Base,  I 
published  June  1978. 


